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2 In this research, inferential theory for hypothesis testing under general 
convex cone alternatives for correlated data is developed. While there exists 
extensive theory for hypothesis testing under smooth cone alternatives with in- 
dependent observations, extension to correlated data under general convex cone 
alternatives remains an open problem. This long-pending problem is addressed 
by (1) establishing that a generalized quasi-score statistic is asymptotically equiv- 
alent to the squared length of the projection of the standard Gaussian vector 
onto the convex cone and (2) showing that the asymptotic null distribution 
of the test statistic is a weighted chi-squared distribution, where the weights 
are mixed volumes of the convex cone and its polar cone. Explicit expressions 
for these weights are derived using the volume-of-tube formula around a con- 
vex manifold in the unit sphere. Furthermore, an asymptotic lower bound is 
constructed for the power of the generalized quasi-score test under a sequence 
of local alternatives in the convex cone. Applications to testing under order 
restricted alternatives for correlated data are illustrated. 

1. Introduction 

Correlated or longitudinal data arise in many areas of science when a response is mea- 
sured at repeated instances on a set of subjects. It is assumed that the measurements on 
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different subjects are independent, while those on individual subjects are correlated with 
an unknown correlation structure (Digglc et al., f 994). In this research, inferential theory 
is developed for hypothesis testing under general convex cone alternatives for correlated 
data using the formula for the volume of a tube around a manifold (curve, surface, etc.) 
on the surface of the unit sphere in an r-dimensional Euclidean space 5R r (Hotelling, 
1939; Weyl, 1939; Naiman, 1990). Testing for order restricted parameters in correlated 
data and testing for a monotone regression become special cases of this general problem. 
Often, interest lies in detecting an order among treatment effects, while simultaneously 
modeling relationships with regression parameters. There exists extensive theory for hy- 
pothesis testing under ordered alternatives with independent observations (Barlow et al., 
1972; Robertson et al., 1988; Silvapulla and Sen, 2004), including smooth cone alterna- 
tives (Takemura and Kuriki, 1997). However, extension of the theory to correlated data 
remains an open problem. 

1.1. Formulation of the testing problem 

Let Yij be the response measured at the jth (j = 1, . . . , n») time point on the ith 
(i = 1,...,N) subject. Let Yj = (Yn, . . . , Yi ni ) T be an n^-dimcnsional vector of re- 
sponse variables. The mean of Y^ is related to the r-dimensional vector of covariates Xy 
corresponding to the r-dimensional parameter vector 7 via a generalized linear model 



where h(-) is the inverse of a link function. We assume that the true distribution is unique 
and all expectations are taken with respect to the true probability measure P. 
The goal is to test the general hypothesis 



where V is an arbitrary finite dimensional vector space of 5R r , r := dim(V © £), £ is a 
closed convex cone with a non-empty interior in W and © denotes the direct or Kronecker 



(1) 




(2) 



Ho : 7 G V against Hi : 7 £ V © £, 7 £ V, 
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sum. Without loss of generality, it is assumed that £ C V x , the orthogonal complement 
of V. Under Tii, dim('y) = r; whereas under Hq, dim(V) < r due to certain constraints 
imposed on the parameters in V. 

Seminal work of Takcmura and Kuriki (1997) has established a solution for the prob- 
lem of testing a simple null hypothesis regarding the multivariate Gaussian mean vector 
A against an arbitrary convex cone alternative for independent observations. In partic- 
ular, they derived the asymptotic null distribution of the likelihood ratio test statistic 
(LRT) for testing 

(3) Hq : A = against Hi : A G JC, 

where K. is a closed convex cone of dimension d with a nonempty interior in ffi (r > d), 
using the techniques of convex analysis. 

1.2. Main results and organization of the article 

The goals of this research include the following. 

1. In Section 2, we derive large-sample properties of the quadratic inference functions, 
extensions of the generalized method of moments (Hansen, 1982), that are required 
for the development of inferential theory with correlated data. 

2. We derive a "generalized quasi-score" (GQS) statistic for the testing problem (2) 
for correlated data in Section 3. Furthermore, it is established that the asymptotic 
null distribution of the GQS statistic, to appropriate statistical order, is equivalent 
to finding the limiting distribution of the squared length of projection of a standard 
Gaussian vector onto the convex cone JC (see Theorem 4) . 

3. In Section 3, we also establish that the asymptotic null distribution of the GQS 
statistic is a weighted chi-squarcd distribution, where the weights are mixed vol- 
umes of tC (Takemura and Kuriki, 1997) and its polar cone (see Theorem 5). De- 
riving computable expressions for the weights in the asymptotic null distribution 
of the test statistic is a tedious and difficult process even for independent data 
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(Takcmura and Kuriki, 1997). Only for special cases, weights are known explicitly 

or can be determined numerically 

4. In Section 4, we express the asymptotic null distribution of the GQS statis- 
tic in terms of certain geometric constants of the volume-of-tubc formula 
(Hotelling, 1939; Weyl, 1939; Adlcr, 1981; Naiman, 1990) for Gaussian ran- 
dom fields (Sicgmund and Worslcy, 1995; Worsley, 1995a,b, 1996) around a convex 
manifold on the surface of the unit sphere (see Theorem 6). We derive explicit 
expressions (in suitable forms for computation) for these geometric constants by 

representing them as integrals over appropriate parts of the manifold. 

5. We derive an asymptotic lower bound for the power of the GQS test under a se- 
quence of local alternatives in K. in Section 5 (see Theorem 9). This lower bound 
demonstrates that the test under restricted alternatives is more powerful than 
the corresponding one under unrestricted alternatives. To the best of the author's 
knowledge, no such lower bound has been derived in the literature even for inde- 
pendent data. 

The article concludes with a discussion in Section 6. 

2. Large-Sample Properties of the Inference Functions 

In this section, we present the large-sample properties of the inference functions 
that are required for later theoretical development of our general testing problem. See 
Pilla and Loader (2005a) for technical details and derivations of the results presented in 
this section. 

Hansen (1982) proposed the generalized method of moments (GMMs) for estimating 
the vector of regression parameters (3 £ B from a set of score functions, where the dimen- 
sion of the score function exceeds that of (3. He established that, under certain regularity 
conditions, the GMM estimator is consistent, asymptotically Gaussian, and asymptoti- 
cally efficient. Qu et al. (2000) extended the GMMs to create a clever approach called the 
"quadratic inference function" (QIF) that implicitly estimates the underlying correlation 
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structure for the analysis of longitudinal data. Their main idea was to assume that the 
inverse of the working correlation matrix, denoted by R _1 (a), is a linear combination 
of several pre-specified basis matrices. That is, R _1 (o:) = X^=i a i Mj, where a%, . . . , a s 
are unknown constants, Mi is the identity matrix of an appropriate dimension and M; 
(I = 2, . . . , s) are pre-specified symmetric matrices with elements taking either or 1 for 
the commonly employed working correlation structures such as exchangeable, AR-1 etc. 

2.1. Properties of extended score functions 

For mathematical exposition, we assume that each subject is observed at a common 
set of times j = 1, . . . , n. Let hj = [/i(X^7), . . . , fo(X?^ 7)] T , where h(X.T is the 
inverse of a link function and the operator V denotes partial derivative with respect to 
the elements of 7; therefore, V hj is the (n x r) matrix (d hi/dji, ... , dhi/dj r ) for each 
i = l,...,N. 

The coefficients a±, . . . , a s in R _1 (q;) are treated as nuisance parameters to create 
the set of subject-specific basic score functions as 



(4) g, (7) := 



Vhf A" 1 / 2 Mi A.71/2 {Yi - hi ) 



vhf a; 1/2 m s a; 1/2 (Y,-h,) 



for i = l,...,N, 



where Aj is the diagonal matrix of marginal covariance of Yj for the ith subject. 

Define the vector of extended score functions for all subjects as gjy(7) := 
Y^iLi Si(7) ■ Note that dim[g Ar (7)] = rs > r = dim(7). The extended score 
vector gjv(7) satisfies the mean zero assumption E7 [15^(7)] = 0, where the expectation 
operator is with respect to the true but unknown distribution of the response matrix Y. 
These estimating equations can be combined optimally using the GMM (Hansen, 1982). 

Let S-y (7) be the true covariance matrix of gi(7), an s-dimensional vector of ex- 
tended score functions defined in (4). We require the following design assumptions for 
deriving the asymptotic theory. 
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Al: The pairs (Y,,Xf) for i = 1,...,N, where X, = (Xji, . . . , Xi„) are (r x n)- 
dimcnsional matrices, are an independent sample from an [n x (r + l)]-dimensional dis- 
tribution F. 

A2: The number of measurements n, on the ith subject is fixed at Uj = n for all 
i = l,...,N. 

A3: The (rs x rs)-dimensional covariance matrix (■y ) := E~^ [gi(7o) Siilo)] ls 
strictly positive definite. 

The independence part of the assumption Al is between different subjects (i.e., with re- 
spect to the index i). The elements of X^ need not be independent of each other; hence, 
this assumption incorporates both time-dependent as well as time-independent covari- 
ates. Moreover, there exists a dependence of on X^ through the link function given 
in (1). 

All throughout this article, E'y denotes an expectation operator with respect 
to the true parameter vector 7 and all expectations are assumed to be finite. Let 
Civil) be the estimator of the second moment matrix of g 1 ('y) so that Cjv('y) := 
A -1 Si=i sfil)- If the mean zero assumption for gjy(7) holds, then A^ 1 Cat(7) 

estimates the covariance oig N ( / y). 

2.2. Fundamental results for the quadratic inference functions 
The quadratic inference function (QIF) is defined as 

(5) Q N (l) :=Ag?;(7)C w 1 ( 7 )gjv(7)- 

If rank of Cat(7) is less than rs or is singular, then the inverse does not exist. However, 
any vector in the null space of Cat (7) must be orthogonal to each of the subject-specific 
score functions gj(7) (i = 1,...,N) and consequently to gjy(7)- Therefore, one can 
replace C^ 1 (7) by any generalized inverse such as the Moore-Penrose generalized inverse. 
Our covariance estimator Cat (7) differs from that of Qu et al. (2000), who define a 
covariance CV with a factor of A -2 , and correspondingly omit the factor of A from the 
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QIF in (5). There are some technical flaws in their results and hence we proceed carefully 
without relying on those asymptotic results. Pilla and Loader (2005a) established precise 
large-sample results for the QIF along with rigorous proofs. The QIF in (5) is minimized 
under unrestricted and restricted spaces to yield the estimators 7 = arg min Qjv(t), 7 = 

7 6 5f 

arg min 2^(7) and 7 = arg min Qn(j), respectively. These estimators can be found 
7ev©<£ jev 

using the iterative reweighted generalized least squares (IRGLS) algorithm developed 
by Loader and Pilla (2006). The IRGLS algorithm avoids the complexity of computing 
the second derivative matrix of QnIi) required for employing the Newton-Raphson 
algorithm recommended by Qu et al. (2000). 

The proof of the next lemma essentially follows from p. 26 of Lee (1996) and hence is 
omitted. 

Lemma 1. Under assumptions A1-A3, the QIF estimator 7 = arg min Qn^) exists 

7e iff 

uniquely and is strongly consistent. That is, 7 — 7 as N — > 00. 

We require the following regularity conditions for further theoretical development. 

A4: The parameter space of 7 denoted by Q C W is compact. 

A5: The parameter space of 7 is identifiable: E'y o [g 1 (7)] 7^ if 7 7^ 7 . 

A6: The true covariance matrix X<y (7) is a continuous function of 7. 

A7: The expectation E^Jgjy^)] exists, finite for all 7 S Q and continuous in 7. 

A8: The subject-specific score functions g^(7) (i = 1, . . . , N) have uniformly continuous 
second-order partial derivatives with respect to the elements of 7. 

A9: The first-order partial derivatives of ^(7) and Cn(j) have finite means and 
variances. 

The importance of assumption A4 is that it enables us to invoke Theorem 1 of Rubin 
(1956); therefore, convergence statements in this article are uniform for 7 in bounded 
sets. 
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Lemma 2. Under assumptions A1-A6, CnIj) — E'y (7 ) as N — > oo. 

Proof. Under assumptions A1-A3 and by the strong law of large numbers, Cat (7) 
converges to its expected value, a non-degenerate limit, for a fixed 7. That is, 

(6) Cjv(7) Efei(7) gf(7)] as A 00. 

Under the stated regularity conditions and Theorem 1 of Rubin (1956), uniform con- 
vergence holds if the compactness assumption A4 holds. Hence, the claim (6) holds uni- 
formly in 7. Lemma 1 combined with the continuity of the function Cat (7) ensures that 
Cat(7) -^-> E[g 1 (7 ) Siilo)} as N — * 00. The claim then follows from the definition of 

S 7o (7o)- ■ 
Let D(7) := E'y[Vg 1 (7)], where Vgi(7) = dg 1 ( / y)/d-y. From the strong law of large 
numbers, Vg Ar (7) — — > D(7) &s N — > 00. This relation combined with Lemma 2, 
enable us to obtain the asymptotic covariance matrix of 7. When there is no ambiguity, 
we drop the subscript 7 and write S _1 (7 ) for the true covariance matrix of gi(7) 
evaluated at 7 . 

Let the estimated covariance of 7 be defined as 



(7) c5v(7) := - [Vg A ;(7)C w 1 (7)Vg Ar (7) 



as N — > 00. 



Lemma 3. Under assumptions A1-A9, N cov(j) [D T (7 ) S^ 1 (7 ) D(7 )] =: 
J _1 (7 ) as N — > 00. 

The proof of the above lemma follows from the previous results. An immediate conse- 
quence of Theorems 3.1 and 3.2 of Hansen (1982) is the next result. The notation 
denotes convergence in distribution. 

Theorem 1. Under assumptions A1-A9, \fN (7 - 7 ) N r [0, J _1 (7o)] as N ^ 

00. 
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Theorem 2. Under assumptions A1-A9, (2^/iV)~ 1 VQA^(7o) 7V r [0,J(7 )] as 
N -> oo. 

2.3. Testing under order restricted alternatives for correlated data 

In the context of correlated data, comparing several treatments, groups or populations 
with respect to their means, medians or location parameters often arise in many areas of 
scientific applications. For instance, one assumes that certain treatments are not worse 
than another. 

The problem of testing under order restricted or constrained hypothesis in longitudinal 
data becomes a special case of (2). Let Y^t be the measurement taken at the jth (j = 
1, . . . , Hit) time point on the ith [i = 1, . . . , nt) subject in the tth (t = 1, . . . , m) treatment 
group. Let N = For mathematical exposition, we assume that ria = n for all 

(i,t) pairs. The mean of Yijt is related to the p-dimensional vector of covariates Xy t , 
corresponding to the p-dimcnsional parameter vector f3 t for the tth group, via E(Yjj t ) := 
h (Xy 4 (3 t + /i t ) , where h(-) is the inverse of a link function, /i t is the treatment effect for 
the tth group. Hence, 7 = (/x T , (3 T ) T with fi = . . . , // TO ) T and (3 — (/3f , . . . , /3^J T . 

The order restricted hypothesis testing problem that is of interest can be formulated 
as Hq : (J, e V against H° : £ <£q, (j, £ Vo, where Vo = {fJ- : fi\ = ■ ■ ■ = fJ, m } and 
€-0 = {fi : Hi > ■ ■ ■ > (J- m } is a particular convex cone. It is clear that r = m(p + 1), d = 
(m — 1) and Vo is the origin of the convex cone Co! hence Vo C Co- This testing problem 
is treated in considerable detail by Pilla et al. (2006). 

3. Hypothesis Testing Under Convex Cone Alternatives for Correlated 
Data 

In this section, we first derive a statistic for the general testing problem (2) using the 
decomposition of 7 € 3? r and next define a new co-ordinate system to transform the null 



10 RAMANI S. PILLA 

space V. Lastly, we derive the asymptotic distribution of the our test statistic under the 
model hypothesis. 

3.1. Generalized quasi-score statistic for correlated data and canonical formulation of 
the testing problem 

Define the generalized quasi-score (GQS) statistic as 

(8) S w := Q N ft) - Q N (j) 

for testing the hypothesis (2), where 7 and 7 are defined in Section 2.2. 

It is more convenient to define a co-ordinate system to transform the null space V. If 
an appropriate transformation is found, we can reduce the general problem to a stan- 
dardized form involving projections of independently and identically distributed standard 
Gaussian random variables as described in the next section. 

Let P be a basis matrix for the space whose columns correspond to the constraints 
imposed by H.q. For example, if al constraints are imposed on Q under 7io, then P is an 
(r x d)-dimcnsional matrix. The choice of the matrix P is problem dependent as shown 
next. 

Lemma 4. The hypothesis testing problem (2) can equivalently be represented in 
terms of the canonical space as testing for 

(9) H% : P T 7 = against H? : P T f £ £j := P T £. 

PROPOSITION 1. Under assumptions A1-A8 and when Hq : P T 7 = holds, 
y/N (P T 7-P T 7) •wZ* ~ N d [0,n{j )} as iV^oo, where ft( 7o ) := P T J- 1 (7 )P. 

Example 1. (Order-restricted testing with three treatment groups). In the case of 
order-restricted testing with three groups, 7 = (/Xi, fj,2, 1^3, fti, P2, Ps) T and V consists of 
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vectors of the form (//, /z, /z, j3\, fo, fa) T which has dimension 4. A basis matrix for V is 



1 -1 



1 1 -2 0, 



so that 



(10) 



P T 7 = 



Mi - M 2 
Mi + M 2 - 2 M 3 



Therefore, the rows of the matrix P span the space V ± . Fig. 1 demonstrates that (i) 
if P T 7 lies in the interior of the convex cone <£i, then P T 7 = P T 7 and (ii) if P T; y 
lies outside <t±, then P T 7 is a projection of P T 7 onto <£% (the orthogonal projection if 

n(7o) = W- 



i 





/ / Cone * 

y p'y 



= P^Y 




Fig. 1. Projection ofP T: y onto the cone ffi /or d = 2, t/ie number of constraints imposed under Hq ■ 



Remark 1. Every 7 £ Q C 3? r admits a unique orthogonal decomposition of the form 
7 = 7i + 7 2 such that 7j £ V and 7 2 <G V and £ C V . It is clear that P T 7 = P T 7 2 ; 
therefore, P T 7 G P T £ since P T V = under H 2 . 

Let Pic Z be the projection of Z onto JC and let || • || denote the vector norm. 
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3.2. Asymptotic equivalence between Sjv and \\Vk Z|| 2 

Owing to Lemma 4, the hypothesis testing problem (2) is equivalent to that of (9). It will 
be established that finding the limiting distribution of Sn for the model hypothesis, to 
appropriate statistical order, is equivalent to finding the limiting distribution of a length 
of a certain projection onto K,. 

We present the two main theorems of this section. The first is based on the quadratic 
approximation to the inference function Qn{i) in an ^-neighborhood of 7 and the 
second is based on the transformed null space. 

Theorem 3. Under assumptions A1-A8, Qn(h)-Qn(i) = N (7 — 7) T J(7o) (7— 
7) + op(l) as N — ► 00. 

PROOF. Let £ := V~N (7 — 70) so that 7 = (7 + N^ 1 / 2 £). From the quadratic ap- 
proximation of the QIF [Theorem 5 of Pilla and Loader (2005a)] the following result 
holds 

(11) Qn (70 + N- 1 ' 2 i) = Q N ( lo ) + 2 (£, W N ) + e J( 7o ) £ + 0P (1), 

where £ <E W is a fixed vector, (, ) is the inner product and Wjv = (2\/~N)~ 1 VQjv(7o)- 
Equivalently, 

(12) Q N (l) - Q N {l ) = 2 (£, Wjv) + f J(7o) I + o P (l). 

The minimizer of the quadratic approximation in (11) is given by = — J _1 (7 )Wjv. 
From Theorem 2, Wjv has a limiting distribution and hence it follows that lies in 
the ball of radius rAr with probability converging to 1. This fact, combined with the 
uniformity of the error term in (11) yields the next result. If £ is the minimizer of the 
quadratic approximation in (11), then the QIF estimator becomes 7 = (7 + iV -1 / 2 £). 
Equivalently, £ = — J _1 (7 ) Wyv + op(l). Therefore, it follows that 

(13) Q N (j) - Q N (7o) = 2 (?, W W ) + f J( 7o ) I + o P (l). 
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Since Wat = — J(7 ) £ + op(l), equations (12) and (13) simplify to Qn(i) — Qn{j) = 

(€ - I) T J(7o) (£-£) + op(1) = N ( 7 - 7 ) T J(7o) (7 - 7) + <*(!)■ ■ 
In the next theorem, we establish the relation between the GQS statistic for corre- 
lated data and the squared length of projection of the standard Gaussian vector Z onto 
K, for independent data. Consequently, we can derive a result for the asymptotic null 
distribution of Sat by employing the seminal work of Takemura and Kuriki (1997). 

Theorem 4. Under assumptions A1-A9, the GQS statistic Sn for testing the hy- 
pothesis under a general convex cone alternative (2) is asymptotically equivalent to the 
statistic \\Vk. Z|| 2 for testing the hypothesis (9), where Z ~ iV r (0, I r ). That is, 

(14) Sat -w \\TkZ\\ 2 as N -> oo. 

Proof. The QIF estimators obtained by minimizing Qn(i) under the spaces V, (V©<£) 
and Q C W are ordered as Qn{^) > 2at(7) > Qjv(7)- From Theorem 3, in an N~ x / 2 - 
ncighborhood of the true parameter vector 7 , Qn(j) — Qn(7) = N (7 — 7) T J(7o) (7 — 
7) + op (I). From Proposition 1, it suffices to consider the transformed hypothesis (9). 
There exists an (r x r)-dimensional matrix L such that L T L = J(7 ); for example, 
Cholcsky factorization of J(7 )- 

It is more convenient to consider <?-paramctrization under the transformed null space 
LV =: V*, where V is the null space under the 7-parametrization. Let — y/Nh'y so 
that = y/NLj, 0= y/NLj and = y/NLj. By Theorem 1, 

(15) (8 - 60) ~» Z ~ N r (0, I r ) as N -> 00. 
Furthermore, Theorem 3 yields 

(16) QNh) = QN{i) + \\8-0\\ 2 + o P {l). 

A given G 3t r admits an orthogonal decomposition as = 0\ + 02 such that 0\ G V* 
and 02 G (V*) . However, orthogonality is not preserved by L; hence (V*) ^ 1 LV ± . 
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The hypothesis (9) can be re-expressed as 

(17) Ho ■ 0i G V*, 2 = against U\ : 1 G V*, 2 G /C and 9 2 0, 

where /C = [L (£ © V)} n {(V*) -1 ] is a d-dimcnsional cone in 5R r since dim^V*)^] = d. 

Similarly, the estimator 6 — ^/NLj admits an orthogonal decomposition as 0i + 02 
such that 0i G V* and 62 G (V*)^. By orthogonality, equation (16) can be expressed as 

(18) Q N [N- 1 / 2 L- 1 (e 1 + e 2 )\ -Q J v(7) = l|0i-0i|| 2 + l|02 - 02|| 2 + op(1). 
Under 7Yq , (18) simplifies to 



(19) 



N- 1 ' 2 L 1 (0! + 2 )] - Qn{i) = ||0i - 0i|| 2 + ||0 2 || 2 + o P (l). 



Under H 3 , 0i G V* and 2 = 0, whereas under Wj 3 , 0i G V* and 2 G K,. 

At first, consider minimizing over Hq : 02 = 0. The right-hand side of (19) is uniquely 
minimized at Q\ = d\ and 02 = 0. By definition, the left hand side is minimized at 
61 =0i. By uniformity of the error term and uniqueness of the minimum, it follows that 
01 = 01 + op(l). Under Tif , the right-hand side of (18) is minimized at Q\ = B\ and 



e-e 




Fig. 2. Projections of estimators of onto the cone K. Asymptotically, the relation = Vic 8 + op(l) 
holds. 



02 = Vk. 02- Therefore, the left hand side of (18) is minimized at 0\ = 0i+op(l) and 02 = 
Vk. 02 + op{l). In effect, minimizing Qjv(-) in (18) under Hq and Hi yield respectively, 
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Q N (l) - Qn(i) = \\02\\ 2 + o P (l) and Q N {j) - Q N {j) = \\0 2 - 2 \\ 2 + o P (l) since 
p! - 0i || 2 = o P (l) and - 0i|| 2 = op(1). Consequently, S N = [Q N (j) - Q N (j)] = 
j|0 2 || 2 -\\02- 9 2 \\ 2 + o P (l) = \\d 2 \\ 2 - \\Vk02 - 02\\ 2 + op(l) = \\VKe 2 \\ 2 + o P (l) since 
[Vk 02 ~02j and Vk 02 arc orthogonal as shown in Fig. 2. Furthermore, since Vk &i = 
under H 3 given that e V* and K. C (V*)" 1 , it follows that S N = \\Vk0\\ 2 + o P (l). 
Moreover, V K O = yielding V K = V K (0- do)- The result follows from (15). ■ 

3.3. Asymptotic null distribution of the generalized quasi-score test 

We derive the asymptotic distribution of the GQS statistic Sn when TC^ holds. Until 
now, we established that (1) the testing problems (2) and (9) are equivalent (Lemma 4, 
Section 3.1) and (2) there exists an asymptotic relation between the statistic Sjv based 
on correlated data for the testing problem (9) and the statistic \\Vk. Z|| 2 based on inde- 
pendent data for the testing problem (3) (Theorem 4, Section 3.2). These main results, 
in conjunction with Theorem 2.1 of Takemura and Kuriki (1997), yield the asymptotic 
null distribution of Sat for the testing problem (9). The weights of this asymptotic null 
distribution are mixed volumes of K. and its polar or dual cone K. (Webster, 1994). 

Let S^ -1 ) be the closed (d — l)-dimensional unit sphere in 5R r , M = K, n be 
the (d — l)-dimensional convex manifold and M° = KP n5' d_1 '. Let tfd-k,k(M, M°) be 
the mixed volumes of M. and Ai° for fe = 0, . . . , d. 

Theorem 5. Under assumptions A1-A9 and when H.q : 7 € V holds, the asymptotic 
distribution of the GQS statistic S n for any c > is 



bedded in 5R r and x 2 is the chi- squared distribution with k degrees of freedom. The \\ f or 
k = is simply a point mass at the origin so that P(xq < c) = 1 for all c > 0. 



lim P(Sat < c) =P(||7>k:Z|| 2 < c) 



(20) 




where d := dim(JC), ld s 



fc-i 



2 7r fc / 2 /r(fc/2) is the volume of the unit sphere S^ k ^ em- 
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The right-hand side of (20) is a weighted mean of several tail probabilities of chi- 
squared distributions; hence, it is often referred to as chi-bar- squared distribution and 
denoted by x 2 (Shapiro, 1988). In general, it is very difficult to derive explicit expressions 
for the weights for the asymptotic null distribution in (20). In certain special cases of JC, 
weights are known explicitly or can be evaluated numerically For instance, for polyhedral 
cones (i.e., the cones defined by a finite number of linear constraints) one can calculate the 
weights. For the general case of non-polyhedral cones, Section 3.5 of Silvapulla and Sen 
(2004) provide a simulation-based approach. 

Takemura and Kuriki (1997) clarify the geometric meaning of the weights when the 
boundary of the cone is smooth or piecewise smooth. However, as the following example 
demonstrates, in the case of order restricted testing problem, one does not have a smooth 
cone or a smooth manifold and hence a more general approach to calculate the weights 
is warranted which is derived in Section 4. 

Example 2. (Asymptotic null distribution of Sat for three treatment groups). Consider 
the problem of order restricted testing with three treatment groups. Under Hq : /i 1 = 
/i 2 = fj, 3 , the asymptotic distribution in (20) has an explicit expression. The convex cone 
JC is the region between the two vectors defining the constraints in (10). Let <f> be the 
angle of the cone JC C 5J 2 at the vertex. Derivation of <f> will be presented in Example 4 
in Section 4.1. Note that the angles (in radians) of JC and JC° at their vertices sum to it. 
Divide the plane into the following four regions: 

1. The cone JC such that Z € JC with a probability of <fi/(2ir) and conditional on 
Z e JC, it follows that VfcZ = Z with \\V K Z|| 2 ~ xl- 

2. The dual cone JC° such that Z = for all Z e JC Q with a probability of [1/2 - 
</>/(2tt)] and conditional on Z e JC°, it follows that \\Vk Z|j 2 ~ xl- 

3. The two regions K) and JC* , where Vic Z is a multiple of one of the two vectors 
defining the constraints in (10); conditional on Z £ JC^ or Z <G /C* with a total 
probability of 1/2, it follows that \\V K Z|| 2 ~ xl ■ 
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For d = 2, the asymptotic distribution in (20) simplifies to a % 2 distribution: 

(21) S N ^\\V K n 2 -{\-^)xl + \xl + ^xl as iV^co. 



Remark 2. The right hand side expression of (21) also occurs in the context of the 
asymptotic null distribution of the LRT statistic for testing for two-component mixture 
model (Lindsay, 1995; Lin, 1997; Pilla and Loader, 2005b). 

Example 3. (Asymptotic null distribution of Sjv under order restricted alternatives). We 
consider the order restricted testing problem discussed in Section 2.3. For convenience, 
we reorder the elements of the parameter vector as 7 = (/x T ,/3 T ) T , with corresponding 
permutations of the rows and columns of J("7 ) [ sce Lemma 3 for the definition of J(7 )]. 
We partition J(7 ) as 

(- 



J(7o) 



V J /3/x J /3/3, 

Let J PL 1 be the appropriate submatrix of J~ 1 (7 ). From the formula for an inverse of 
a partitioned matrix, it follows that J A* A* — ^J/i/x — J fi{5^ flfi^ (3fi) ' Since subjects 
in different groups are independent, the variance matrix J^t 1 is diagonal. Let Z> ~ 
N m (0, Q), where Q is a pre-specified diagonal variance matrix. 

In the order restricted testing problem, the asymptotic null distribution has an ex- 
plicit expression. Under H§ , for any c > 0, the result (20) reduces to lim P (Sn < c) = 
J2k=o p{d—k + l, d\ Q) P (Xd-k - c )' w h ere Q = andp(fc, d; Q) is the level probabil- 
ity that the projection of onto £0 with a weight vector Q consists of exactly k distinct 
points [Section 2.4, Robertson ct al. (1988)]. The unknown weight vector is replaced with 
Q. This problem is developed and treated in considerable detail by Pilla et al. (2006). 
The weights p(k,d; Q) are also referred to as x 2 weights (Robertson ct al., 1988). Such 
weights appear in the null asymptotic or exact distribution of several test statistics when 
there are inequality constraints on parameters. In certain cases, exact expressions for 
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these weights are available and in other cases, one may obtain approximations or bounds 
(Silvapulla and Sen, 2004). 



4. Asymptotic Null Distribution of S^: The Volume-of-Tube Formula 

In this section, we derive explicit expressions for the weights in the asymptotic null 
distribution of Sat in (20) by representing M. and /C in parametric form and in turn using 
the Hotelling-Weyl-Naiman volume-of-tube formula. 

4.1. Parametric representation of M and K, 

From a geometrical perspective, the rows of the (r x d)-dimensional matrix P span the 
space V . From the orthogonal decomposition in Section 3, it follows that 7 = 7! +72 = 
7 j + Pu such that j 1 £ V and for some u £ 3? d . Let TV := {u : Pu £ €}. Given that 
7 £ <t = {Pu: u £ Af}, the hypothesis (2) can be expressed as 



In order to represent IC and M. in a parametric form, we require the following result. 

PROPOSITION 2. The matrix P* := (L _1 ) T P forms a basis for (V*)" 1 , where V* = 
LV. 

Proof. Let y £ V* so that y = Lz for some z £ V. Let x £ C(P*), the column space of 
the matrix P*, so that x £ (L _1 ) T Pu for some vector u = (u 1 , . . . , u d ) T £ W l . It follows 
that, (x,y) = (u T P T L- 1 )(Lz) = (Pu, z) = since z £ V and Pu £ V^. Therefore, the 
space spanned by P* is a subspace of (V*) ± . However, this subspace and (V*) ± have the 
same dimension d, therefore they must be equal. ■ 
We return to the orthogonal decomposition of 9 = 6\ + 62 such that 6\ £ V* and 
62 £ (V*) x . Proposition 2 ensures the following representation: 62 = 7-Vy*)-L 6 = V-p* = 



(22) 



Hq : u = against : u £ TV. 



P* [(P*) T P*1 (P*) T 9 = P* [P T J(7 )P] 1 P 7 ^-^ since L T L = J 



( 7o ). There- 
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fore, 



2 = VNP* [P t J( 7o )P] 1 P T (7i+Pu) 



(23) 



VNUu since P T J 1 = 0, 



where 



(24) 



H := P* [P T J( 7o ) P]" 1 P T P. 



Next, the d-dimensional cone K, = {L (<£ V)} n {(V*)^} and the (d — l)-dimensional 
manifold M. = /CnS^ 1 ' are represented in the parametric form by considering a vector 
function T(u) :Afc3t d ->Mc K^" 1 ), where ueAf. That is, 



where H is an (r x d) matrix defined in (24). In effect, we can redefine AA ;= {T(u) € 

: u £ A/" C and /C := {H u: u E JV C 3? d }. 
Example 4. (Explicit Expressions for the convex cone Af C 3fJ 2 and the angle </> of /C). 
We return to the problem of testing under order restricted hypothesis for correlated data. 
If 7i° : > n 2 > fj, 3 , then the choice of P is given in Example 1. Therefore, -y 2 = P u 
corresponds to /ij ~ (u + u 1 + u 2 ),/i 2 = (u — u 1 + u 2 ) and /i 3 = (u — 2it 2 ), where 
u = +/i 2 +/i 3 )/3. The constraints fi 1 > fi 2 and fi 2 > fi a under Tl° yield respectively, 
u 1 > and u, > u 1 /3. Consequently, the convex cone Af = {u: u 1 > 0, u, > u 1 /3}. It 
is clear that u lies in the cone bounded by the vectors vi = (1, 1/3) T and V2 = (0, 1) T . 
The cone JC is then bounded by the 02-component of LPvi and L P Hence 



yields an explicit expression for the angle <fi defined in Example 2. 

4.2. Asymptotic null distribution of Sat in terms of the geometry of A4 

As a first step, we establish the connection between the distribution of a squared length 
of projection of Z onto K, and the volume-of-tube problem. We take a different approach 



(25) 



T(u) :={Hu:ueAAc3fi d and||Hu|| = 1}, 
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from Lin and Lindsay (1997) in order to cast the problem in the general framework of 
this article. 

The geodesic (or angular) distance between two points on any manifold is defined 
as the shortest measured distance between the points within the manifold itself. Let 
X(g, M.) or X(</>, M.) be the spherical tube around the topological (d — l)-dimensional 
manifold M. of Euclidean radius g or geodesic radius <fi embedded in iS( d_1 ), where 
g = \/2[l — cos(0)], u G N and dim(A/") = d. Since iS>( d_1 ' is also a manifold and 
X C 6> ( - <i ~ 1 " ) , the geodesic distance between two points on X is the length of the segment 
of the great circle (arc) connecting the two points. We view each ray {(rj : ( > 0} as 
a cone on which to make a projection, yielding Zjy that depends on rj. We redefine the 
cone as JC :— {(rj: £ > 0, |jr7|| = 1} to yield Vic Z = arg min \\(rj — Z|| 2 = sup (77, Z) + , 
where (■, ■) , denotes the positive part of the inner product. In effect, we have 



(26) 



\VkZ\ 



sup (T7,Z), 
rjeic 



In order to reduce the problem to that of a uniform process, we condition on ||Z|| 2 
and integrate over the conditional distribution. Consequently, from (26) we can express 



P(||^Z|| 2 <c) = 



SUP { T] 



< 







( sup ( V ,U) + ) < - 


l|Z|| 2 = z 



fr(z) dz 



(27) 



^sup^(j7,U) + < \J- ] f r (z)dz, 



where U = (Zi/||Z||, . . . , Z,./||Z||) is uniformly distributed on S r , an r-dimcnsional unit 
sphere embedded in $t'- r+1 \ and f r (z) is a x 2 density with r degrees of freedom. Therefore, 
the right-hand side of (20) can be determined from (27), provided the probability in the 
integrand can be found, at least approximately. 

The uniformity property of U reduces the problem of finding P(sup^ (77, U) + < \J c/z) 
to that of determining the volume of the tube 1(g,A4) including the end points correc- 
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tions proposed by Naiman (1990). Consequently. 



(28) 



sup (rj,U) + < J- 



) 



P[U G T(</>,.M)] 



where w r _ 1 = 2 7r r / 2 /r(r/2) is the volume of iS r and $x(</>) is the volume of 1(g, ./Vf). 
Therefore, (27) and (28) establish a connection between \[Pk, Z|| 2 and volume of the tube 
1((/>,A4) around M. embedded in S r . Essentially, we established that the distribution of 
jlP/cZH 2 can be determined explicitly by finding #m(0) which equals cos~ 1 (y/c/z) for 
any < 4> < tt/2. 

The asymptotic expansion of the tail probability of the sup (?7, Z) + can also be ob- 
tained using the Euler-P oincarc characteristic £ method, developed by Adler (1981) 
and Worslcy (1995a,b), where the expectation of the £ of an excursion set is evaluated. 
Takcmura and Kuriki (2002) establish the equivalence between the tube and Euler char- 
acteristic methods under the assumption that A4 is a manifold with piecewise smooth 
boundary. 

We motivate the geometric concepts through the order restricted alternatives for cor- 
related data. We define a corner to mean a point where two faces of the boundary of A4 
meet. 

Example 5. (Geometry of M. for the order restricted alternatives). We assume that 
rii = n (i = 1, . . . , N) and the number of subjects in each group is equal so that we have a 
balanced design. In this case, lo = 2 and u) 1 = 2 7T. First consider three treatment groups 
(i.e., m = 3), then the number of restrictions d equals two corresponding to fi ± < /i 2 < fx 3 . 
Therefore, dim(/C) = 2 and M. is just an arc with two end points. Suppose m — 4 corre- 
sponding to three constraints, then A4 is a spherical triangle. The interior corresponds 
to /itj < fj, 2 < /Lt 3 < /a 4 with three corners fi 1 = fx 2 = /i 3 , /Lt 2 = [i 3 = (i i , \i x = fj, 2 and 
/Lt 3 = /i 4 and three edges \i 1 = /i 2 , /i 2 = /i 3 , ^ 3 = /.t 4 . 

Example 5 demonstrates that determining , &m{4') f° r d>3 depends on the geometry 
of M. Naiman (1990) derived expressions for the volume of a tube by decomposing the 
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tube into different sections, corresponding to the main part of the manifold, hemispher- 
ical caps along boundaries of the manifold, circular wedges at the boundaries and so 
on. Adding up these terms yields a scries involving partial beta functions [Lemma 3.6 
of Naiman (1990)]. Substituting these terms into (27) yields a series involving partial 
gamma functions; the first four terms of which are given in the next theorem whose 
proof essentially follows from Pilla and Loader (2005b) and hence is omitted. 

Theorem 6. Under assumptions A1-A9 and when Ha : J E V holds, the asymptotic 
distribution o/Sat for a d-dimensional manifold A4 and for any c > is given by 

lim P(Sjv > c) =P(\\V K Z\\ 2 > c) 

N—>oo 

(29) + (4 + Vl + T) P {xl-a >c)+o (c( d - 5 )/ 2 e- c / 2 ) as c -> oo, 

where k is the (d — 1)- dimensional volume of the manifold M, k 2 is the measure of 
curvature of M, £q is the {d — 2) -dimensional volume of the boundaries of M, l\ is the 
measure of rotation of the boundary, £2 is the measure of curvature similar to k 2 , v is 
the measure of rotation angles at points ( or along edges ) where two boundary faces meet, 
v 1 is the combination of these rotation angles with the rotation of the edges, and r is the 
measure of the size of wedges at corners where three boundary faces of A4 meet. 

If the manifold Ai is a single point (i.e., d = 1), the result (29) simplifies 
to 0.5P(xi>c). If M is one-dimensional (i.e., d = 2), the result reduces to 
(« /27r)P(x! > c) + (4/4)P(xf > c), where k q is the length of M and £ is the 
number of boundary caps which equals 2. This last result is same as that obtained by 
Lin and Lindsay (1997); however, we provide an explicit formula for k based on the 
parametric representation of M. which is derived in the next section. 
Remark 3. For convex manifold M., n 2 = in the asymptotic expansion (29). In the 
case of order restricted alternatives, the manifold M is a high-dimensional tetrahedron 
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whose corners correspond to the constraints imposed on 7 under Hi; hence v ^ 0. Also, 
£2 = in the case of order restricted testing problem. 

Owing to Theorem 6, the weights in Theorem 2.1 of Takemura and Kuriki (1997) 
for the independent data case can be determined using the volumc-of-tubc formula. 
When the critical geodesic radius of the manifold (Naiman, 1990) is greater than or 
equal to 7r/2, all of the coefficients in (29) are nonnegative leading to a finite mixture 
of chi-squared distribution or x 2 distribution. The critical radius is greater than equal 
to 7r/2 if and only if the smallest cone containing the manifold is convex. Lemma 2.1 of 
Takemura and Kuriki (2002) provides a formula for computing the critical radius. 

4.3. Explicit expressions for the geometric constants 
We derive explicit expressions (in suitable forms for computation) for the geometric 
constants in (29) using the representation of T(u), defined in (25), and its derivatives. 
Our main goal is to reduce the evaluation of the constants to integrals over appropriate 
parts of M. 

The profound result of Gauss-Bonnet theorem (Do Carmo, 1976; Milman and Parker, 
1977) connecting curvatures of manifolds with the Euler-Poincarc characteristic 
(Worsley, 1995a,b; Adler and Taylor, 2004) can be employed to find some of the 
geometric constants appearing in Theorem 6. When A4 is two-dimensional (i.e., d = 3), 
the number of pieces contributing to A4 minus the number of holes equals £ . In particu- 
lar, k 2 + 1\ + v = 2 7r£ — k which eliminates the need to compute k 2 , £\ and vq directly. 
Remark 4- Lin and Lindsay (1997) assume that the cone is convex and smooth; hence 
no corners (i.e., v = 0). Both their Theorem 3.1 and the result for d = 3 presented 
in Section 4 of Lin and Lindsay (1997) become special cases of our general result 
established in Theorem 6. 

In the parametric representation of A4, the function T(u) has an embedded con- 
straint ||Huj| = 1 for a d-dimensional vector u £ Af. This embedded constraint means 
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that the manifold M. is of dimension (d — 1). Hence, we can represent in terms of a 
(d — l)-dimensional parameter vector p with u = u(p) and express T(p) = T[u(p)]. For 
example, such a transformation can be carried out using the polar co-ordinates. Denote 
the domain of p by A/"*. We express T(p) = [T 1 (p 1 , . . . , p d _ 1 ), ■ ■ ■ ,T r (p 1 , . . . , p d _ 1 )] for 
the parametric representation of M . 

A10: The transformation T(p) is one-to-one and each T l (I = l,...,r) is twice 
continuously differentiable on Af* C ^RS d ~ 1 '. 



The following expressions are derived under A10. Define an [r x (d— 1)] matrix S(p) 
[Tr(p) • ■•T d _ 1 (p)], where T fc (p) = 9T(p)/9p fc for fc = 1, . . . , (d - 1) and Ti(p), . 
T,j_i(p) are the column vectors. 

The volume of the manifold .M is expressed as 



«o = / dct[S T (p)S(p)] 1/2 dp. 
Joe W* 



IpeM 

Finding £2 is essentially similar to that of finding K2 by simply treating each face of the 
boundary as a new manifold. Therefore, we describe the method to find k 2 , although it 
is zero for convex manifolds. In the case of testing under order restricted alternatives, 
the constant £2 = 0. The measure of curvature of M. is expressed as 

K a = [ \ [T(p) - (d - l)(d - 2)] dct [S T (p) S(p)] 1/2 dp, 
JpeAf* 1 

where T(p) =T£%Z^ Ya=1 {™kk CT « ~ m kl with 

d 



zuf k = ef [S T (p) S(p) 



dp 



■S(P) 



I d _i-S(p) {S T (p)S(p)}^ S T (p) 



and e; as the basis vector for 5ft^ -1 ^. 

The volume of dA4, the boundary of .M, is £q while measures the curvature 
of d A4; both of these need to be determined for each face of the boundary. For in- 
stance, on the face where p d _ 1 = 1, define Sf(p) := [Ti(p) • • • Td_2(p)] and B<j_i(p) := 



Id— 1 — S(p) {S T (p) S(p)} 1 S T (p) Td_i(p), where <p is a normalizing constant. 
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and 



to= det[Sj(p)S t (p)Y , Up 
JpedM 



f £r(p) dct[S[(p)S t (p)] 1/2 dp, 
J pedM 



where 



d-2 



Mp) = -£el [Sf(p) st(p)]- 1 



k=l 



sT(p) 



B d _i(p). 



Similarly, 



/ u,(p) det[Sf(p)S t (p)] 1/2 dp, 

J oed 2 M 



iped 2 M 

where d 2 M. is the region or corner at which two boundary faces of M. meet, 

d-3 



= [sf(p) st(p)]- 1 



fc=l 



Sf(p) 



[B d _ 2 (p) + B d _i(p)] tan 



and 4>{p) is the angle between 13^-2 (p) and Bd-i(p). Since .M is of dimension (d — 1), 
d M. and 9 2 M. arc of dimensions (d — 2) and (d — 3), respectively. Lastly, we consider 
the edges where two boundary faces meet. If we consider the edge where p d _ 2 = p d _ 1 = 1 
and define Sj(p) := [Ti(p) ■ • ■ T d _ 3 (p)], then 

u = / v B (p) det[Sf(p)S t (p)] V2 dp, 

Jpea 2 M 

where v (p) = arc cos[(B d _ 2 (p), B d _i(p))]. 

In order to find r, we need to calculate the area of the spherical triangle which is 
achieved by the Euler's formula: area of the triangle equals (<f>i +4>2+4>3~ 7T ), where 0i, 2 
and 03 are the three internal angles of the triangle. Loader and Pilla (2007) describe the 
method of determining these angles by first finding the vectors defining the corners of 
the triangles. 
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5. Power Under a Sequence of Local Alternatives 

In this section, we derive an asymptotic lower bound for the power of the GQS statistic 
under a sequence of local alternatives in /C. This plays an important role in comparing 
the result with a test against the unrestricted alternative. To the best of the author's 
knowledge, a lower bound has not been established in the literature even for independent 
data; hence it would be an interesting one to derive. 

From the parameterization defined in Section 4.1, we can express 7 = ~f 1 + P u such 
that 7j € V and u g AT. Following the hypothesis (22), we consider a sequence of local 
alternatives of the form 

u* 

(30) u„ = -= for a fixed vector u* g TV C 5R d and N = l,2,.... 

From the derivation of 9 2 in (23), the relation O2 = ^fN~tiu N = Hu* holds under the 
sequence of local alternatives (30) . 

As a first step, we define a statistic for testing 7io :76V against the unrestricted 
alternative H2 '■ 7 g Q as Sjy := [Qn^j) — Qn{'j)), where 7 and 7 are defined in Section 
2.2. Using the arguments similar to Theorem 1 and the result (15), one can establish 
that 7} ~ 7Y r (Hu*, l r ) under the sequence of local alternatives. 

The arguments given in Robertson et al. (1988) and Pilla and Loader (2005a) yield 
the following result. 

Theorem 7. The asymptotic local power of the unrestricted test statistic S^r for a 
sequence of alternatives (30) is 

(31) lim P(S^>6 1 )=P[x^(^ 2 )>6 1 ], 

N—>OQ 

where b\ > is a constant, S = ||Hu*|| and x^(<5 2 ) is the chi-square distribution with a 
non-centrality parameter S 2 and with r degrees of freedom. 
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The above result is equivalent to Theorem 7 of Pilla and Loader (2005a); however, 
here the non-centrality parameter S is represented in terms of H. Theorem 7 yields an 
exact local power and the next one establishes a lower bound for S^. 

Theorem 8. A lower bound for the asymptotic power of S^, under a sequence of 
alternatives defined in (30), is lim P (S^ > 61) > 1 — $ {\fb~x — 5) . 

It is worth noting that finding the asymptotic power for S n under the sequence of local 
alternatives in K, is hard and it does not have a simple weighted non-central chi-squared 
distribution with a pre-specified non-centrality parameter. The following result gives an 
asymptotic lower bound for the power of Sjv under a sequence of local alternatives (30). 
It demonstrates that Sat under restricted alternatives (i.e., testing for H.q against Ti.\) is 
locally more powerful than S* N under no restriction (i.e., testing for TIq against 0.2)- 

Theorem 9. A lower bound for the asymptotic power o/Sat for a sequence of alter- 
natives (30) is lim P (Sat > b 2 ) > P (N(6, 1) > Vh) = 1 - $ (Vfo - S), where b 2 > is 
a constant, S := j|Hu*|| and $ is the standard Gaussian cumulative distribution. 

PROOF. Let d = arg min HZ* — b|| 2 , where Z* ~ f r (Hu*,I r ). Consequently, d = Vk, Z*. 
be/c 

Further let, C := {cHu* : c> 0} and d* = arg min ||Z* - a|| 2 , then d* = V c Z*. Since 

as £ 

£ C K, it follows that ||d|| > ||d*||. Equivalcntly, \\V K > \\Vc Z* || - 

Let a = Hu*/||Hu*|| so that ||a|| = 1. It is clear that j|ef || = (a, Z*) The proof of 
(a,Z*) ~ N(S, 1) is presented next. We have E [(a, Z*)] = (Hu*,Hu*) /||Hu*|| = 6 and 
variance V(||d*||) = j|H u*|| ~ 2 V[(Hu*) T Zt] = 1, since Z* ~ iV r (Hu*,I r ). Therefore, 
the relation (a, Z*) ~ N(S,1) holds. Furthermore, ||'?- > /c Z* || > (a, Z*), since the right- 
hand side is the length of the projection of Z* onto H u* G K, and (a, Z*) is distributed 
as N(5, 1). The theorem claim follows from (14). ■ 
Example 6. (Comparison of the local power of and Sat for the order restricted 
testing problem). Define ~f N := (fi 1 n ,/j, 2 N , . . . , /Lt m N , f3±, . . . , f3 m ) T as a sequence of 
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local parameter vectors for N = 1,2, . . .. The local alternatives take the form [i k N = 
[i h _ x N + e k _ 1 /^/N for k = 2, . . . ,m, where t 1 , . . . , e m _ 1 arc fixed negative constants. As 
N — > oo, the sequence of local alternatives approach the null hypothesis Tio : 7 € V. 

Consider three treatment groups (i.e., m = 3) leading to six possible orderings, with 
each order corresponding to an arc on the unit circle. Union of these six arcs comprises the 
unit circle S . Due to the balanced design assumption, each of these arcs is of the same 
length; therefore, the angle of the cone /Cis</> = 7r/3 = 60°. At the level of significance 
a = 0.05, the critical values corresponding to the two tests Sat and are b 2 = 3.820 
and bi = 5.991, respectively. Tabic 1 presents the asymptotic lower bounds on the local 
power for the two tests. The table also presents the exact asymptotic local power for 
S^r obtained using the asymptotic formula (31). It is clear that except for S = 0, the 
asymptotic local power of Sat is better than that of S^, in terms of both the lower bound 
and the exact power. 

Table 1 

Comparison of the local power of the tests for a given 8, the non-centrality parameter, when 
m = 3 and <j> = it/ 3 = 60°. 



Test 


6 

1 2 3 4 5 


Sat Lower bound 

Exact 
S^r Lower bound 


0.025 0.170 0.518 0.852 0.980 0.999 
0.050 0.133 0.416 0.771 0.957 0.996 
0.007 0.074 0.327 0.710 0.940 0.995 



6. Discussion 

In this research, inferential theory is developed for the problem of testing under convex 
cone alternatives for correlated data. Such a problem occurs when interest lies in detecting 
ordering of treatment effects, while simultaneously modeling relationships with other 
covariates. The testing problem (2) is also applicable to the analysis of clustered multi- 
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categorical data. In this framework, Ya = (lilt, . . . 1 Y%Kt) T denotes the if-categorical 
response on the ith observation in the cluster t, where = 1 if category j (J = 1, . . . , K) 
is observed and otherwise. 

We established that the GQS statistic is asymptotically equivalent to the squared 
length of the projection of the standard Gaussian vector onto an arbitrary convex cone 
with a nonempty interior. We further derived the asymptotic null distribution of the GQS 
statistic under convex cone alternatives for correlated data as a weighted chi-squarcd dis- 
tribution. The weights in the asymptotic distribution are the mixed volumes of the convex 
cone and its polar cone which do not have explicit expressions except in special cases. 
For non-polyhedral cones, closed-form expressions for the weights are very complicated 
and therefore; often a simulation approach is employed for computing them [Section 3.5, 
Silvapulla and Sen (2004)]. In this article, explicit formulas are derived for the calculation 
of these weights using the Hotclling-Wcyl-Naiman volume-of-tube formula. 

Furthermore, an asymptotic lower bound is derived for the power of the test under 
a sequence of local alternatives in K. for correlated data which establishes that the test 
under restricted alternatives is more powerful than the test under no restriction. Note 
that Barlow et al. (1972) and Robertson et al. (1988) derive the asymptotic power only 
under specified alternative hypothesis. The current theory is applicable to many practical 
problems of interest including testing for a monotone regression function and for the 
analysis of clustered multi-categorical data. 
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